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Abstract 

We study the high speed collision and reconnection of Abrikosov-Nielsen-Olesen cosmic strings in the type-II regime of 
the Abelian Higgs model, that is, scalar-to-gauge mass ratios larger than one. Qualitatively new phenomena such as multiple 
reconnections and clustering of small scale structure have been observed in the deep type-II regime and reported in a previous 
paper, as well as the fact that the previously observed "loop" that mediates the second intercommutation is only a loop for 
sufficiently large mass ratios. Here we give a more detailed account of our study, which involves 3D numerical simulations with 
the parameter f3 = m 2 scalar / rn 2 gauge in the range 1 < f3 < 64, the largest value simulated to date, as well as 2D simulations of 
vortex-antivortex head-on collisions to understand their possible relation to the new 3D phenomena. 

Our simulations give further support to the idea that Abelian Higgs strings never pass through each other, even at ultrarela- 
tivistic speeds, unless this is the result of a double reconnection; and that the critical velocity for double reconnection goes down 
with increasing mass ratio, but energy conservation suggests a lower bound around 0.77 c. We discuss the qualitative change in 
the intermediate state observed for large mass ratios. We relate it to a similar change in the outcome of 2D vortex-antivortex 
collisions in the form of radiating bound states, whereas we find no evidence of the back-to-back reemergence reported in 
previous studies. 

In the deep type-II regime the angular dependence of the critical speed for double reconnection does not seem to conform to 
the semi-analytic predictions based on the Nambu-Goto approximation. We can model the high angle collisions reasonably well 
by incorporating the effect of core interactions, and the torque they produce on the approaching strings, into the Nambu-Goto 
description of the collision. An interesting, counterintuitive aspect is that the effective collision angle is smaller (not larger) as 
a result of the torque. Our results suggest differences in network evolution and radiation output with respect to the predictions 
based on Nambu-Goto or j3 — 1 Abelian Higgs dynamics. 



Introduction 

The discovery of cosmic strings, first proposed by Kib- 
ble pQ would revolutionize our understanding of particle 
physics at the extremely high energies present in the very 
early Universe. It could signal a 'superconducting' phase 
transition and give information on the particle interac- 
tions before the transition, or it might provide the first 
evidence of superstring theory. So far there is no evi- 
dence of their existence, but strings can lead to a wealth 
of detectable astrophysical phenomena and there is an 
increasing number of surveys and searches looking for 
observable signatures. These include Cosmic Microwave 
Background (CMB) anisotropies, gravitational lensing, 
wakes, gravitational radiation (GW), cosmic rays and 
gamma ray bursts, among others. Gravitational effects 
are determined by the adimensional parameter Gfi which 
for most models of cosmological interest falls in the range 
1CU 13 — 1CP 6 (we use units with h = c = 1 throughout; 
G is Newton's constant and /i the mass per unit length 
of the strings). Strings with higher mass per unit length 
are already ruled out by these observations, see the clas- 
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sic reviews PHS] and the more recent updates [BtiTT] , and 
references therein. 

The formation of cosmic strings and supcrstrings is 
a generic outcome of cosmological phase transitions [1 
and of some inflationary models, in particular those 
based on Grand Unified theories [12l [13] . More recently 
it has been appreciated that some models of brane 
inflation could also lead to cosmic superstrings Q31 [T5] . 
Once formed, a string network is expected to reach a 
scaling solution in which statistical properties such as 
the distance between strings or the persistence length 
become a fixed fraction of the horizon size (the age of the 
Universe). The energy density in strings decreases with 
time but the expansion of the Universe pumps energy 
into the string network -by increasing the contribution 
of long strings- and this is balanced by energy losses to 
radiation. If these are efficient enough, the contribution 
from the strings to the energy density of the Universe 
remains a small, constant fraction of the dominant 
form of energy (matter or radiation) and is potentially 
observable. Radiation is emitted by oscillating loops, 
formed when a string self-intersects, and in bursts. The 
latter are produced by cusps (sections of the string 
which acquire near-luminal speeds), by the final stages 
of collapsing loops and, to a lesser extent, by kinks 
created when strings reconnect. The reconnection, or 
intercommutation of strings is therefore an essential 
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process that determines and maintains the long-term 
scaling behaviour of the string network. 

In this paper, a companion to |16j . we focus on the 
high speed intercommutation of Abelian Higgs strings in 
the deep type-II regime. The terminology is borrowed 
from superconductors, where type-I (type-II) indicates 
a critical parameter /3 < 1 ((3 > 1); (3 is the ratio of 
scalar to gauge excitation masses, squared (see below). 
An important difference between type-I and type-II is 
the interaction energy between parallel vortices, which 
is attractive for type-I and repulsive for type-II [17]. 
Cosmic string intercommutation was first investigated 
numerically by Shellard for global strings |18) and 
later by Matzner for type-II (3 — 2 Abelian Higgs 
strings [19]. They pointed out that in ultra-relativistic 
collisions there is a critical center-of-mass velocity v c 
(depending on the collision angle) beyond which strings 
pass through: a loop forms, expanding rapidly from the 
collision point, that catches up with the reconnected 
strings and produces a second intercommutation. A 
more recent study |20j focused on double reconnection 
and showed that it proceeds differently in type-I and 
type-II strings, with the loop only forming in type-II 
collisions. In fact, one of the main results in [T5] that 
we elaborate on here is that the loop only forms in 
deep type-II (/3 >> 1) collisions. Although we do not 
determine the precise value of (3 for which the transition 
occurs, it is somewhere between j3 = 8 and j3 = 16. For 
/3 > 16 we see the loop forming, while type-II collisions 
with j3 < 8 produce instead a blob of radiation that can 
look like a loop but is not (see later sections and [16]). 
Another important point is the angular dependence of 
v c as a function of the collision angle. In ref. [2U] it was 
shown that it is dictated by the geometry and speeds 
of the strings after the first reconnection, which can be 
calculated in the Nambu-Goto approximation. In the 
deep type-II regime we expect core interactions to play 
an increasingly important role, even before the collision, 
and this is indeed what we will report here. 

The simulations in [20 had decreasing resolution with 
increasing (3, and only explored two values of (3 in the 
type-II regime {(3 = 8, 32), but the results suggested that 
the critical velocity for the second reconnection would go 
down as a function of (3. This is interesting because one 
of the distinguishing features of cosmic superstrings, as 
opposed to the Abrikosov-Nielsen-Olesen (ANO) strings 
[2T1 |2"2"] considered here, is their low intercommutation 
probability P ~ 1CT 3 - 10" 1 23J, which leads to different 
scaling properties, in particular to denser networks (al- 
though in an expanding background the effect is weaker 
than the pt 2 ~ \jP dependence one might expect for 
the density p at cosmic time t |24|). However, if the 
critical velocity of strongly type-II Abelian Higgs strings 
decreases to the extent that it becomes comparable to 
the average velocity of the network, the effective inter- 
commutation probability could be much less than one. 



This was one of the motivations behind the present study. 

Our results confirm the claim of [50] that the critical 
velocity for the second reconnection goes down as a 
function of (3, although energy conservation suggests 
this decrease cannot go on indefinitely. We will return to 
this point later. Furthermore, while studying the critical 
velocity, we found multiple intercommutations. That is, 
processes where the strings exchange ends three times or 
more. We only found multiple reconnections for /3 > 16. 
We will show that this is related to a qualitative change 
in the nature of the intermediate state, from localized 
radiation to a string loop, which determines the process 
leading to the second intercommutation. For 1 < f3 < 8 
we found that the previously reported loop is just an 
expanding blob of radiation, while for (3 > 16 we find 
a topological loop rapidly expanding from the collision 
point. 

An interesting question that was not addressed in |16] 
is whether these multiple reconnections and the blob- 
to-loop transition could be related to, for instance, two- 
dimensional bound states of the vortex-antivortex system 
or whether they are a purely three-dimensional effect. In 
order to shed light on this connection we also simulated 
the high speed head-on collision of a vortex and an an- 
tivortex in 2D. Abelian Higgs vortex-antivortex scatter- 
ing in two dimensions was studied years ago by Myers, 
Rebbi and Strilka[25] (see also [26]). They reported that, 
beyond a certain critical speed, the vortices reemerge and 
the way in which they reemerge depends on (3: for f3 < 4 
they bounce back and for (3 > 8 they pass through (and 
for velocities below the critical speed the pair annihi- 
lates into radiation). For large (3 we agree with their 
results: we found that, for (3 > 6.4, the vortices reemerge 
as if they have passed through. However, for (3 < 6.2 we 
find no evidence of backscatter; instead, we always find a 
bound radiating state. The energy profile looks as if the 
vortex-antivortex have passed through but a closer look 
at the magnetic field shows they are just fluctuations. 
They cannot escape each other's influence and radiation 
keeps being emitted from the area around the collision 
point until all the energy has been radiated away. As we 
will discuss, these two behaviours are consistent with the 
change we observe in the intermediate state in 3D col- 
lisions. The "pass through" behaviour at high (3 corre- 
sponds to a 3D loop while the 2D radiating bound states 
at low (3 correspond to the 3D "blob". On the other 
hand, the relation with multiple reconnections is more 
subtle, and it appears that the number of reconnections 
is mainly a three-dimensional effect. 



A few comments on cosmic string intercommuta- 
tion 

String intercommutation poses some interesting 
puzzles from a theoretical point of view. ANO strings 
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are topological, they cannot break. When two strings 
segments collide there are, in principle, three possible 
outcomes: 

a) They can exchange partners and reconnect (in- 
tercommute). The reconnected strings have a slightly 
shorter length. This is the default outcome and, as we 
shall discuss, seems to occur always when Abelian Higgs 
strings meet. 

b) In near parallel collisions of Type-I strings at low 
velocity the attractive interaction makes them stick 
together, e.g. type-I Abelian Higgs strings[371 125] . 
The network then has junctions between strings with 
different winding numbers (and different energy per 
unit length). String networks with junctions have been 
extensively studied analytically in the Nambu-Goto 
approximation [2"9l 130] and in this regime they give an 
extremely good fit to numerical simulations |3IJ. 

c) The strings can simply pass through each other. 
The relative probability of outcomes a) and c) is 
an important distinguishing feature of ANO strings 
versus cosmic (super)strings. The intercommutation 
of fundamental superstrings, and of D-strings, is a 
quantum process and, as such, has some probability 
of not happening[32 . This is also an expected feature 
of higher dimensional models -in which the apparent 
collision of the strings can be a four dimensional illusion, 
the strings are actually not intersecting in the higher 
dimensional space- or models with extra internal de- 
grees of freedom [5] [33] . Cosmic superstrings have been 
shown to give reconnection probabilities (P) as low as 

10-3 _ io-![23]. 

An interesting observation is that one would naively 
expect outcome c) to be the result of any sufficiently fast 
collision of ANO strings, since the natural timescales 
of the microscopic fields can be much slower than the 
string crossing time 34 . In the case of domain walls, 
for instance, the scalar field profiles pass through each 
other undistorted in ultrarelativistic collisions 35 . The 
head-on collision of a vortex and an antivortex in 2D 
at ultrarelativistic speeds is also known to result in 
the two passing through each other [25] • However this 
"free passage" is not observed in 3D string collisions. 
In the Abelian Higgs model, all (numerical) evidence to 
date pSl HH [2Qj points to the conclusion that P = 1: 
Abrikosov-Nielsen-Olesen (ANO) strings with unit 
winding always reconnect at least once, even at 
ultrarelativistic speeds (even in case b) there is some 
evidence that the strings will first reconnect and then 
settle into a junction [3T]). What is observed, instead, 
is that beyond a certain critical collision velocity v c , 
the strings may reconnect a second time and effectively 
pass through with some distorsion. In terms of network 
evolution this is as if the strings have not reconnected, 
and one can talk of an effective intercommutation prob- 
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FIG. 1: From ref. |20| (a) Initial positions and orientations of 
the strings in the center of mass frame. The strings lie in the 
x = const, planes and approach each other with speed v. The 
arrows indicate the orientations of the strings, which form an 
angle a. (b) The configurations after one intercommutation. 
If v ~ c, the kinks' motion along the strings is negligible and 
the connecting horizontal segments are practically antiparallel 
and immobile, making a second interaction possible. 



ability P e ff being less than I. An important question 
is to model the dependence of the critical velocity for 
non-reconnection as a function of collision speed and 
angle, and to understand its dependence with /3. 

At lower speeds, an interesting class of analytic results 
on intercommutation is based on the moduli space 
approximation [36 . In the Bogomolnyi limit (/3 = 1) 
[57] . the moduli space or geodesic approximation gives 
a very good description of slow vortex-vortex collisions 
in two dimensions and predicts the observed right-angle 
scattering[3S]. Although these arguments are only valid 
at low speeds and in the absence of core interactions, 
they lead to the expectation that in three dimensional 
collisions reconnection is inevitable. Consider the 
configuration in fig. [T^), two strings colliding with angle 
a and velocity v. If the incoming strings lie in the yz 
plane and move in the x-direction, the collision looks in 
the xz plane like two vortices that will scatter at 90° 
[2"6l [551 [3"9"] while in the xy plane it is a vortex-antivortex 
collision, which we expect will annihilate. In 3D the 
energy from this annihilation goes into creating the 
connecting segments and the rest is emitted as radiation. 

This idea has been extended to other models such as 
D-strings [ID] and non- Abelian strings [Ml EI])- In many 
cases intercommutation is expected with probability one 
in the regime of validity of the moduli space approx- 
imation, which -as we just emphasized- requires two 
conditions: low speed, so higher order frequency modes 
are not excited, and negligible core interactions (such as 
in near parallel collisions not far from the Bogomolnyi 
limit) . 

It is worth stressing that even if these conditions 
are satisfied before the collision, they will not hold 
after intercommutation because the "new" portions of 
string that are generated between the receding strings 
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are necessarily almost antiparallel around the point 
of collision, and their core interactions are crucial in 
understanding what happens next. This is probably 
the reason why a prediction for the angular dependence 
of the critical speed in [3D] fails to agree with the 
data from numerical simulation. This prediction is 
based on an energy argument which does not take into 
account the interaction between the string cores, and 
this approximation fails immediately after the collision. 

By using the thin string approximation and taking into 
account the effect of core interactions after the collision, 
ref. [SDJ obtained a semi-analytic expression for the criti- 
cal velocity for double reconnection v c as a function of the 
collision velocity v and angle a. This works quite well for 
type-I and moderately type-II ANO strings. Consider the 
situation in figure [T]b) after the first intercommutation. 
Within the Nambu-Goto approximation -thus, ignoring 
core interactions- we can express the angle 6 between 
the horizontal segments and the speed w at which the 
horizontal segments move apart in terms of v and a: 



w = sin(a/2)/7(w) 



cos(<ty2) 



cos(a/2)/(vj(v)) 
y/1 + (cos(a/2)/( U7 ( V )))2 



(1) 



(2) 



Notice that for high collision speeds the horizontal 
segments are almost antiparallel, S ~ n and they move 
slowly w ~ 0. If they are not receding very fast, the 
anti-aligned segments will be attracted to each other 
and will annihilate, causing the strings to reconnect a 
second time. In the unphysical limit of a collision at the 
speed of light, the bridging segments would not move 
at all and would be antiparallel, so the second recon- 
nection is expected with probability one. The angular 
dependence of the threshold speed v c is then found with 
equations ([I]) and |2| if one assumes that reconnection 
will happen below a threshold "escape" speed w t and 
above a threshold angle S t , close to antiparallel. The 
threshold values wt and 8t are left as free parameters 
and determined by the best fit to the numerical data. 
While this model works well near the Bogomolnyi 
limit P ~ 1, we will show that it does not work so 
well in the deep type-II regime studied here, indicat- 
ing that core interactions before the collision are also 
important. We will return to this point in the discussion. 

Equation shows that fast collisions will produce 
very slowly moving connecting segments. Conversely, 
a slow collision usually creates a highly curved region 
after intercommutation which will accelerate under its 
own tension and acquire large speeds. This is one of 
the mechanisms that helps the string network maintain 
a typical speed < v 2 > ~ 0.5, even in an expanding 
Universe. 



Simulations 

The Abelian Higgs model is the relativistic version of 
the Ginzburg-Landau model of superconductivity. It is 
described by the lagrangian 

(3) 

where <f> is a complex scalar field and A^ is a U(l) 
gauge field with field strength F^ — d^A v — d v A^ 
(fi, v = 0,1, 2, 3). 

The ground state has |</>| = r\ and zero electric 
and magnetic field. The fluctuations about this 
vacuum define two mass scales: the scalar exci- 
tations have m sca i ar — yXr] and the gauge field 
excitations, m gauge = \[2er\. Classically, the only 
relevant parameter in the dynamics is their ratio, 
P = (m scalar /m gauge ) 2 = A/2e 2 , which also character- 
izes the internal structure of the ANO vortices. 

Magnetic cores repel and scalar cores attract, so the 
interaction between vortices is determined by which of 
these cores is the largest: parallel ANO vortices repel 
for P > 1 and attract for P < 1. The Bogomolnyi limit 
P = 1 is a critical value where both effects cancel and 
parallel vortices do not interact. In this paper we are 
interested in the P > 1 regime, analogous to a type-II 
superconductor, and in this case the vortices have an 
inner "scalar" core of radius ~ TO 7cLar m which the 
scalar field departs from its vacuum value and vanishes 
at the center. This is surrounded by a larger, "gauge" 



core of radius 



gauge 



where the magnetic field is 



non-zero. The repulsive interaction produces a torque 
that tends to anti-align two colliding strings. This will 
play an important role later. 

Here, as in jTD], we follow the numerical strategy 
of [IH1 [2D]: we use a lattice discretization and place a 
superposition of two oppositely moving ANO strings on 
a three dimensional lattice. This configuration is evolved 
using a leapfrog algorithm. The initial configuration 
is determined by two parameters: the center-of-mass 
speed v of the strings when they are far apart and the 
angle a between them (every collision can be brought to 
this form by an appropriate Lorentz transformation [18] ) . 

We also impose "freely moving" boundary conditions: 
after each round the fields inside the box are updated 
using the equations of motion, and the fields on the 
boundaries are calculated assuming the strings move 
unperturbed and at constant speeds at the boundaries. 

All 3D simulations were done on a 200 x 200 x 400 
grid. Unless otherwise stated, we use a lattice spacing 
a = 0.2 and time steps At — 0.02, so the Courant 
condition (here At < a/yS) holds. Some simulations 
were repeated with a = 0.1, in particular those with 
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P = 4, to confirm the results. 

Our simulations are optimized for the deep type-II 
regime. By solving the two-dimensional, static vortex 
equations one finds that in a static straight cosmic string 
about half of the potential energy in the scalar core is 
contained within a radius V / 2w~ c ^ ;al ./(/3), where / is 
a slowly varying function with /(l) = l,/(64) = 1.4. 
Lorentz-contraction gives an extra factor 7(f)" 1 in the 
direction of approach, with j(v) = 1/vl — v 2 . This 
is the smallest length scale that has to be resolved. 
Without loss of generality we take A = 2, 77 = 1, which 
ties the unit of length (and time) to m~ calar = l/v2- 
The scalar core is resolved by at least three lattice 
points up to a center of mass speed of v ps 0.94 — 0.96, 
which is indicated explicitly in the diagrams in figure 
[2] (it is higher in the /3 = 4 simulations because those 
have a = 0.1). The initial string separation is fixed to 
^>V^P7( v ) m 7 C alar = 5\/]3j(v). This would be about 
five times the actual core radii for /3 = 1 but as j3 
increases, the core sizes increase and for large j3 one has 
to check that the gauge cores do not overlap in the initial 
configuration. For j3 — 64 the overlap in total energy 
from the tails of the gauge cores is less than 1% when 
7(1;) = 1, that is, when calculated on static vortices. 
Finally, the two-dimensional simulations of vortex- 
antivortex head-on collision and reemergence were done 
on a 800x800 grid with a lattice spacing of 0.1 and 
At = 0.02. In this way we resolve the vortex cores with 
at least 3 points up to a speed of v ps 0.985. The initial 
separation is the same as in the 3D simulations. We 
used absorbing boundary conditions. In all simulations 
(2D and 3D) energy is conserved to better than 5% until 
the radiation hits the boundary (which determines the 
dynamical range). 



Results 

We simulated the collision of cosmic strings at 
13 = 1,3.9,4.0,4.1,8,16,31,32,33,49 and 64 for various 
speeds v and angles a to find the threshold velocity 
above which the strings effectively pass through each 
other. The results for selected values of /3 are shown in 
figure [2] Results for f3 = 3.9 and 4.1 were qualitatively 
similar to /3 = 4; also, results for j3 — 31, 32 and 33 were 
qualitatively similar. Some salient features were already 
reported in reference |16j . 

The first thing that is apparent in fig. [2] is that the 
minimum critical velocity for a second reconnection goes 
down as a function of /3, in agreement with what was 
observed in [2U]. The dependence of this lowest critical 
velocity v c>m i n with /3 and the range of collision angles 
for which it is observed is seen in table U 

Secondly, there is a new phenomenon of multiple 
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4 


0.92 


42° 


8 


0.92 


62° 


16 


0.90 


90° - 110° 


32 


0.88 


80° - 120° 


49 


0.88 


85° - 125° 


64 


0.86 


105° - 120° 



TABLE I: The lowest value of the critical velocity for double 
reconnection as a function of f3, and the collision angles at 
which it is observed. The range in a indicates the possible 
existence of a plateau in the critical velocity. 



intercommutations, which appears to be related to a 
change in the nature of the intermediate state from a 
non-topological blob of radiation to a loop. 

The images of the intercommutation process in figures 
[3j [4] and [6] show isosurfaces of the scalar field with 
\4>\ = 0.4. At this value only about 20% of the potential 
energy is contained within them, but it allows us to 
visualize the evolution of the Higgs field most effectively. 
A tube twice the radius would contain about 60% of the 
energy (to be precise, a threshold of \(j)\ = 0.8, which has 
twice the thickness of the tubes shown, contains 62% of 
the scalar potential energy for j3 = 16, 57% for /3 = 32 
and 52% for £ = 64). 

It is clear from these images that the strings do 
not always intercommute once or twice, as previously 
observed, but also three and four times for particular 
values of the initial speed v and angle a. An odd number 
of reconnections results in overall intercommutation 
of the strings, and an even number in the strings 
effectively passing through, so we can still speak of 
a threshold velocity for the strings passing through. 
However, each reconnection creates small structure on 
the strings in the form of a left- and a right-moving 
kink. In some cases it is not easy to distinguish be- 
tween one and three reconnections, or between two 
and four reconnections, just by looking at the energy 
isosurfaces in the intermediate state. But the resulting 
kinks are clearly visible in the final state and can be 
counted; in case of doubt we use this criterion (see fig. |4| . 

Successive reconnections therefore lead to left- or 
right-moving "kink trains" , groups of up to four closely 
spaced kinks (one for each intercommutation). The 
inter-kink distance within these trains is a few core 
widths, at the time of formation (see figs [3] and |4| . 

A multiple intercommutation process for type II 
strings unfolds as follows. After the collision in which 
the strings exchange ends for the first time two things 
can happen: 
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• For /3 > 16, an expanding loop forms after a short 
delay. If v > v c , the loop catches up with the two reced- 
ing strings and these reconnect again through the loop. 
This creates a highly curved central region in each string 
(sometimes, for the lower collision angles, the loop is not 
as sharp as in figure [3] and the bridge is very pronounced, 
making the intermediate state look more like a junction). 
The central regions will move towards each other and 
in some situations mediate a third reconnection (see fig. 
[3]). After the third reconnection there are two almost 
antialigned string segments, and if they are receding 
sufficiently slowly a fourth reconnection is possible. This 
is the largest number of intercommutations we have 
seen, and only for (3 = 16. On the other hand, if v < v c 
the second reconnection does not take place because the 
string loop does not catch up with the strings. In this 
case it will contract again and eventually decay into 
radiation, sometimes after one or a few oscillations. 

Triple intercommutations are quite generic for collision 
speeds and angles on the boundary between the regions 
in parameter space where we observe one and two 
reconnections. For f3 = 64 we see a few triple intercom- 
mutations in a window around v ~ 0.87, a ~ 94°. The 
box is already somewhat small and we expect a larger box 
with increased dynamical range would show more multi- 
ple intercommutations, but this remains to be confirmed. 

• For 1 < f3 < 8 the energy isosurfaces look somewhat 
similar, but they reveal a very different intermediate 
state. The "loop" in figure [6] is just a blob of ra- 
diation with no topological features: the (covariant) 
phase of the Higgs field around the "vortex" that 
makes the loop shows no net winding around the loop, 
/q" : (dg4> — ieAg(j))d6 = 0. This is clearly visible in 
the third timestep, where the "loop" meets the string 
bridges, breaks and is absorbed - a real loop of string 
would not be able to break if it carried a net winding. 
In this case the maximum number of reconnections is 



two. The blob slows down the receding strings (thereby 
lowering the critical velocity) and facilitates the second 
reconnection. However, whether or not the second 
intercomutation takes place is still determined by the 
string bridge (see fig. [lj. We therefore see, as expected, 
a good agreement between the data and eqs. [Tf2]m fig. 
[2] for P = 4 and 8. 

We now turn to the simulations of the 2D head-on 
collision of a vortex and an antivortex at ultrarelativistic 
speeds. The parameters (/3, v) of the simulations we 
performed are listed in table [TTJ While we confirm the 
general picture of ref. [55], we have slightly different 
results for the state after the collision: 

• For j3 < 6.2, the emerging vortex and antivortex 
settle in a bound, oscillating state, which completely de- 
cays into radiation for all speeds between 0.9 < v < 0.98. 
This behaviour extends to the /3 < 4 regime studied by 
Myers et al. [25] . who interpreted the outcome as the 
reemergence, back-to-back, of the original pair. We see 
no evidence of this reemergence. A typical configuration 
after the collision is shown in figure [5] Although 
some timesteps could be mistaken for back-to-back 
reemergence of the vortex-antivortex pair, subsequent 
evolution makes it clear it is only localized radiation. 

Note also that, since the vortex-antivortex pair never 
reforms after the collision, the expectation of back-to- 
back reemergence, suggested in [H] by analogy with the 
global vortex case, also does not apply. 

• For (3 > 6.4 and high collision speeds the pair 
reemerges as if they passed through, and the critical 
velocity above which the vortex-antivortex pair passes 
through each other goes down with increasing j3, from 
around v = 0.98 for /? = 6.4 to (/3, v): (6.6, 0.98), (6.8, 
0.98), (7, 0.95), (8, 0.95) and (32, < 0.7). This agrees 
with ref. [25]. 
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FIG. 2: The number of intercommutations for a range of collision velocities and angles for j3 — 4, 8, 16, 32, 49 and 64. 
The symbols A, +,• and * stand for 1,2,3 or 4 intercommutations respectively. The dotted lines show a two-parameter fit to 
simulations with collision angles below 150°, based on the Nambu-Goto approximation 20] The dashed line is a one-parameter 
fit adapted to the deep type-II regime at high collision angle (see text for explanation and fit parameters) . Simulations above 
the horizontal line resolve the scalar core size by less than three lattice points and are therefore less reliable. The point at 
(/? = 32, q = 137.6, v = 0.92), indicated by a square is inconclusive because the intercommutation happens just beyond the 
dynamical range. The /3 = 4 simulation has lattice spacing a = 0.1, all others have a — 0.2 
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TABLE II: The parameters (/3,v) of the 2D simulations of 
vortex-antivortex collisions described in the text. The simu- 
lations in boldface are those where the vortex-antivortex pair 
reemerges as if they had passed through. In all other cases 
the pair annihilates into radiation, sometimes after forming 
a short-lived pulsating bound state. In the simulations in 
parentheses, the vav-pair annihilate directly into radiation. 
The case /3 = 4, v = 0.95, indicated by a star, is shown in fig. 




Discussion 



Our results suggest some interesting differences be- 
tween the (3 » 1 regime and the much more studied 
/3 = 1 regime when it comes to the intercommutation be- 
haviour and the resulting small scale structure. Some of 
these differences can be traced back to the core interac- 
tions, in particular the repulsion between parallel, deep 
type-II strings, which in a 3D setting can distort the im- 
pact angles and velocities of the strings. The angular 
distorsions can be parametrized and could in principle 
be incorporated in numerical simulations of cosmic string 
networks and analytic studies of small scale structure. 




FIG. 3: Isosurfaces of the scalar field with \4>\/rj < 0.4 for a 
collision with (/3 = 32, a = 122.7, v — 0.88) showing a triple 
intercommutation. From left to right, up to down: snapshots 
at t — 3, 4.7, 5, 6.9, 10, 17, 20, 22. Time is measured in inverse 
scalar masses (see text). At t = 4.7 the strings collide. Notice 
the distorsion around the point of collision (see also figure FF|. 
After the first intercommutation a loop emerges at t = 6.9. 
The loop catches up with the receding strings and intercom- 
mutes at the connection points (t = 10, second intercommu- 
tation), creating a highly curved central region in each string 
(t = 17). These move towards each other and produce a third 
intercommutation at t = 20. Two sets of three closely spaced, 
left- and right-moving kinks are visible on each of the strings 
in the last panel (indicated by arrows on the upper string). 
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FIG. 4: After the third intercommutation there are two, al- 
most anti-aligned, string segments close to each other (see 
the last panel of figure [3] If these are receding slowly, even 
a fourth intercommutation is possible. This figure shows the 
typical string configuration after the fourth intercommuta- 
tion. The arrows show four closely spaced kinks moving up 
each string; the two corresponding down-moving trains are 
visible below the collision point. (f3 — 16, a = 122.7, v = 
0.92), same isosurfaces. 




Multiple reconnections and the nature of the inter- 
mediate state 



As reported in [16], we have observed a qualitative 
change in the process determining the second intercom- 
mutation. For 1 < (3 < 8 we see the emergence of a blob 
of radiation after the first intercommutation. A blob can- 
not cause a second intercommutation by itself, it is ab- 
sorbed when it reaches the strings. In this case, whether 
or not the second intercommutation takes place is deter- 
mined by the geometry after the first intercommutation. 
This situation is well described by eqs. []][2]as found in 
|20j . On the other hand, for /3 > 16 we see the emergence 
of an expanding topological loop that mediates the sec- 
ond intercommutation. We compare this transition with 
the reemergence of the vortex-antivortex pair in highly 
relativistic two-dimensional collisions: as j3 increases, the 
critical speed for 2D vortex-antivortex reemergence goes 
down below the critical speed for double reconnection in 
3D (and even below the 3D universal bound on the crit- 
ical speed v ~ 0.77 [To] , see next subsection). Therefore, 
for sufficiently large /3, reemergence of the v-av pair is un- 
avoidable. In 3D the reemergent vortex-antivortex pair 
leads to a loop. Thus, for j3 > 16, the outcome is always 
the same, the loop always forms. However, whether or 
not this loop leads to a second reconnection depends on 
whether the loop catches up with the receding strings, 
which is, in turn, determined by the velocity and angle 
of the strings before the collision: for high v there is a 
second reconnection, for lower v there is only one recon- 
nection. It is also in this loop-mediated regime, of which 
the lower bound is between /3 ~ 8 and f3 = 16 according 
to our simulations, where the multiple or higher order in- 
tercommutations take place. By contrast, for low /3, the 
critical speed in 2D for the reemergence of the v-av pair 
is so high that we do not see the emergence of a string 
loop in 3D. Instead we see what we describe as a radi- 
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FIG. 5: A typical 2D head-on collision of a vortex and an- 
tivortex in the regime where where the outcome is a radiat- 
ing bound state (f3 below 6.2). The figure shows the mag- 
netic field of a vortex and antivortex collision with P = 4 
and v = 0.95. All other parameters and timescales are like 
in the previous 3D simulations. The Grey to Blue colorscale 
indicates the magnetic field pointing out of the page, while 
Black to Red indicates the magnetic field points into the 
page. Black (grey) indicate points whose magnetic field is 
> 20% of the maximum B value, which is attained at the 
core of the incoming vortex (antivortex). The intermediate 
colours red/orange/yellow/light green (dark blue/blue/light 
blue/dark green) are in decreasing steps of 5%. From left 
to right and top to bottom the panels show snapshots at t = 
0,7,9,10,12,14,17 and 18. First, we see the vortex and antivor- 
tex at t=0. Second, we see a "pair" reform (t=7), but they 
are actually two blobs of radiation. In the third figure (t=9) 
we see a configuration which could be mistaken for a back- 
to-back reemergence. However, we then see a second blob 
forming (t=10), and the first pair fizzling out. The magnetic 
field oscillates around zero. At t=12 the second "pair" (blob) 
forms and some radiation is exchanged. By this time, the 
magnetic field of the first "pair" has already switched sign. 
This proves that the first "pair" had no topology. In the 
sixth panel (t=14) polarity is opposite to that in the second 





FIG. 6: After the first intercommutation a radiation blob 
emerges. The blob catches up with the bridge and is ab- 
sorbed. However, before it is absorbed, a loop seems to be 
formed. This "loop" is not topological, it breaks resulting 
in "Dracula's teeth", (j) = 4, a = 122.7, v = 0.98). The in- 
teraction between the blob and the strings slows them down, 
facilitating a second reconnection, but -unlike in figure |]-the 
blob by itself cannot mediate this second reconnection. From 
left to right, up to down: snapshots at t = 2, 4.5, 6.5, 8. 
Note the antialignment of the bridging segments between the 
strings, as in figure [I] 



ation blob, the 3D equivalent of the bound, oscillating 
state in 2D. An interesting open question is to identify 
more precisely the value of j3 (between 8 and 16) at which 
the transition between these two regimes occurs in 3D. 



The critical velocity for double reconnection: /3 de- 
pendence, angular dependence and a universal lower 
bound 



The dependence of v c with the collision angle shows 
that, as we go deeper into the type-II regime, core 
interactions are playing an important role, especially at 
high collision angles. With better resolution and more 
data points than in [20] we actually see a difference with 
the f3 w 1 behavior. The best fits with eqs. [T][2] do not 
quite agree with the data for large j3 (e.g. /? > 16). 
These fits underestimate the critical velocity for impact 
angles close to antiparallel (or else do not properly 
account for impact angles below 90°). Discarding data 
points with impact angles larger than 150° leads to the 
fits shown in figure [2] with the dotted line. We found 
for j3 = 4, 8 the following parameters: Wt = 0.202,0.238 
and 5 t = 136.4°, 139.3°. For (3 = 16,32,49 
and 64 we found w t = 0.328,0.357,0.261,0.408, 



S t = 144.0°, 136.1°, 139.0°, 135.7° respectively, (we note 
that for (3 = 49 the fit has fewer relevant data points). 

The deviation between the fits and the data at high 
collision angle is, by itself, not so surprising since 
in the deep type-II regime two things change: First, 
the appearance of the loop means the angle between 
the string segments is not expected to be relevant. 
The second intercommutation will occur if the loop 
touches the receding segments, irrespective of their 
mutual angle. Second, the core interaction affects 
the state before the collision, it produces a torque 
that deforms the strings and tends to anti-align the 
colliding segments. This means the true collision angle is 
actually larger than the initial value. The torque results 
from the attraction between vortex and antivortex 
and vortex-vortex repulsion in the orthogonal plane, 
which compete in the type-I regime but add up in type-II. 

One might expect that a modification of the NG fit 
of [20] to account for this offset in collision angle should 
give a better fit in the deep type-II regime: 



w = sin[(a — cvq)/2] 



1 



cos(<S/2) 



cos[(q - a )/2]/(wy(v)) 
Vl + (cos[(a-a )/2]/M«))) 2 



(4) 



(5) 



with ctQ negative. This expectation is not borne out by 
the data. In fact, if anything, the data prefer a modifi- 
cation with a smaller angle (positive ao - see table |III[ ) . 
And this can be understood by looking at figure [7j Al- 
though the actual collision angle (the angle at the point of 
collision) is larger than a, the motion of the strings after 
the collision is determined by an effective collision angle 
£ that is smaller than a. This is because the antialign- 
ment causes a larger portion of string to be annihilated 
so, after a quick transient, the strings look as if they had 
collided with angle £ < a and the bridge segments (see 
fig[lj are further apart than they would have been in the 
absence of a torque. 

If the collision is sufficiently close to antiparallel, 
the relation between the angle of approach, a, and 
the effective collision angle £ is universal (see fig. [7|): 
tana/2 = (1 + sin£/2)/ cos£/2, that is, £ = 2a — it, or 
ao = 7r — a, expected to be valid for large a. . This 
one-parameter fit is shown with a dashed line in figure 

The fit parameter is w cr i t and it has values (f3,w cr i t ): 
(16, 0.233), (32, 0.243), (49, 0.251) and (64, 0.275). 



The offset in collision angle (see table III) for (3 > 16 



indicates that the torque is strong even at very high colli- 
sion speeds, and suggests strong distortions for low speed 
collisions. This is potentially very interesting from the 
point of view of the radiation coming from cosmic strings. 
Usually, the radiation bursts from reconnection are sub- 
dominant to those from cusps and kinks [121 03]. But 
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64 


0.257 


234.7 


58.7 



strings in flat space, to 0.63 — 0.51 in field theory 
simulations with j3 = 1 with cosmic expansion [45 | 146 ) ). 
So, double intercommutations may be less rare than in 
type-I collisions, but they are still rare events. 



Summary and outlook 



TABLE III: The best fit parameters of eqs. [4] and [5] A pos- 
itive (negative) offset ao indicates the actual collision angle 
is smaller (larger) than in the initial configuration. We see 
that ao becomes positive for /3 > 16 and remains positive and 
relatively large. 




FIG. 7: The bending of the strings at the point of collision 
due to the intervortex potential, modelled in the figure as a 
sinusoidal perturbation, a is the initial collision angle when 
the strings are approaching but still far apart. The torque 
tends to anti-align the segments around the collision point (see 
the second panel of figure [3j. Subsequently, a vertical segment 
of length ~ 2L will annihilate from each string, leaving the 
receding strings as if they had collided with angle £. Note 
that, although the torque tends to increase the collision angle, 
the effective reconnection angle £ is actually smaller than the 
initial collision angle a. 

here reconnection bursts are enhanced because longer 
segments of string are annihilated. The amount of ra- 
diation produced by cosmic string reconnection might 
therefore be somewhat larger in deep type II collisions 
than would be expected from analytic, Nambu-Goto ar- 
guments. 

Finally, our results appear to confirm the claim [2U] 
that the critical speed goes down with increasing j3, from 
v c ~ 0.96 at P = 1 to v c ~ 0.86 at (3 = 64, although this 
reduction cannot go on indefinitely. A crude (universal) 
lower bound for v c follows from energy conservation |16j : 
if the strings anti-align locally and a portion L of each 
string is annihilated, the maximum energy available 
to the loop is 2L"f = 2irR, with R the maximum loop 
radius. If R < L/2 the second intercommutation cannot 
take place, which happens for v < \fl — 4/-7T 2 ~ 0.77. 
These values of v c are to be contrasted with the average 
root mean square velocity of the network, which has 
not yet been investigated for type-II strings (known 
values range from v rms = l/\/2 ~ 0.71 for Nambu-Goto 



We have investigated numerically the intercommu- 
tation of Abelian Higgs strings in the deep type-II 
regime for selected values of j3 = rn 2 scalar /m 2 gauge » 1 
up to (3 — 64, the highest value studied to date. Our 
study shows interesting differences with the (3 = 1 
behaviour and raises some puzzling questions. New 
effects arise due to the strong interactions of the string 
cores. Multiple reconnections were already reported in 
[16) . and also a qualitative change in the intermediate 
state after the first reconnection, with truly topological 
loops only appearing for /3 > 16. It is also in this 
regime (/3 > 16) where we find the higher order (three 
or more) intercommutations. Further work is needed 
to understand if the window closes for (3 > 64. We see 
fewer multiple intercommutations, but this could be 
simply due to the limitations of the dynamical range. 

As j3 increases, we find a lower critical velocity for 
double (or multiple) reconnections, in agreement with 
[2"U] , but for very large (3 we are unable to describe the 
angular dependence in detail. The fit derived in [20] 
works well for j3 < 8 and maybe even for (3 — 16 but 
not in the deep type-II regime (3 > 32. The interaction 
between the magnetic cores produces a torque that 
tends to anti-align the string segments, so the actual 
collision angle is higher than in the initial configuration, 
but this is only true around the collision point, and this 
portion of string quickly disappears. In fact, due to this 
antialignment, the collision results in the annihilation of 
a larger segment of string so in fact the strings behave 
as if they that had collided with a smaller angle. The 
simplest way to model the bending of the string, as a 
sinusoidal excitation, is shown in fig. [7] and gives the 
angle £ in terms of a: £/2 = a — n/2. Using £/2 instead 
of a/2 in eq. [I] we find a one-parameter fit to our data 
with initial collision angles higher than 120 degrees. The 
fit is shown with a dashed line in figure [2] 

Even for lower collision angles, the fit derived in |20) . 
with or without an adjustment to include the torque 
between the colliding strings, does not fully describe 
the angle dependence of the critical velocity. Our data 
(see fig. [2]) is possibly better described in the deep 
type-II regime by a curve of opposite concavity with 
a plateau between collision angles of 80° and 120°. A 
more detailed study is necessary to understand the angle 
dependence of the critical velocity. 

We also simulated two dimensional vortex-antivortex 
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head-on collisions in an attempt to understand the 
new 3D effects. We confirmed the result of Myers 
et al. [35] that for high enough collision speed the 
vortex-antivortex pair reemerges some time after the 
collision. The critical speed for reemergence goes 
down with increasing /3, and for (3 = 32 is already 
lower than v ~ 0.7. On the other hand, for (3 < 8 
the velocity needed for the vortex-antivortex pair to 
reemerge is so high (v > 0.98) that we only see a 
(bound) radiating state. In particular, for f3 < 4 Myers 
et al. report the backscatter of the vortex-antivortex 
pair whereas we always see a radiating bound state. 
This is however not inconsistent with [25] given their 
much smaller dynamical range (see fig. [5]). So we 
see a transition between reemergence in forward direc- 
tion vs a bound radiating state, we never see backscatter. 

Further work is needed to determine the critical 
f3, in 2D head-on v-av collisions, that distinguishes 
the regime where a bound radiating state forms from 
the regime where the v-av pair reemerge as if they 
had passed through. We can locate this transition 
somewhere around 6.2 < f3 < 6.4 but the main problem 
in determing this critical j3 is the high collision speed 
(v > 0.98) needed for reemergence, which leads to very 
bad resolution. 

Our expectation is that the 2D radiating bound states 
at low (3 should roughly correspond to the 3D blob, and 
that 2D forward reemergence should correspond to the 
3D loop. For not too large (3 (in particular below 6.2) 
the critical velocity for passing through in 2D is so high 
(> 0.98c ) that we do not probe it in either the 2D or 
3D simulations, we always see a radiation blob. As f3 
increases, the critical velocity for passing through in 2D 
goes down and at high j3 (in particular >32) the critical 
velocity for forward reemergence is so low (< 0.7c) that 
all the 3D simulations are in this regime, and all show 
loop formation. 

Although suggestive, our interpretation of the results 
is not conclusive. To confirm this picture one should 
study the transition region 8 < (3 < 32, 0.7 < v < 0.95 
both in 2D and 3D and verify the extent of this 
correlation (it is important to note that we have no 
data between j3 = 8 and (3 — 16). Also, we do not 
necessarily expect a precise, one-to-one correspondence 
in the critical values of f3 and v because the energy 
requirements to reform a vortex-antivortex pair in 2D 
are different from those needed to form a loop in 3D. 

The effect on the cosmological signatures of strings 
is hard to predict, as the stronger core interactions in 
the deep type-II regime affect energy loss mechanisms 
in several, competing, ways [16;. In general, we expect 
a relative enhancement of the radiation contribution 
from kinks and reconnections at the expense of cusps 



(suppressed by the kinks [171135]) and loops (suppressed 
by the lower critical velocity for double reconnection) . 

Regarding kinks, a new feature is the presence of kink 
trains resulting from multiple reconnections. These are 
rare, but once a kink train is formed its decay time is 
comparable to that of a single kink, and because of its 
microscopic size (it is only a few core widths in length) 
it is very unlikely to be disrupted by intercommutation 
with another string segment. We conclude that the small 
scale structure of strongly type-II Abelian Higgs string 
networks could be somewhat more clustered than the 
predictions based on the Nambu-Goto approximation 
with P = 1 [4T)H5T] . although nowhere near the prolif- 
eration of kinks expected in a network with iunctions|52j. 

Regarding reconnections, they are enhanced in two 
ways. First, while we confirm that strongly type-II 
ANO strings always reconnect (P — 1), an effective 
intercommutation probability P e // < 1 due to multiple 
reconnections will still lead to denser networks and 
therefore more reconnections. Second, most of these will 
be at low velocity, and we have argued that antialign- 
ment plays a role even in relativistic speed collisions, 
so we expect a strong effect in low velocity collisions. 
So we would expect stronger and more frequent bursts 
of radiation and cosmic rays than for other string 
types (lower (3 and also superstrings) where reconnec- 
tion bursts are always negligible or subdominant [331 133] • 

Further work is needed to understand these effects 
quantitatively, and how they affect the cosmological 
bounds. But the upshot of the work presented here is 
that core interactions are expected to cause significant 
differences with respect to the predictions from both 
Nambu-Goto strings and field theory Abelian Higgs 
strings in the Bogomolnyi limit. 

Finally, our results confirm once again that Abelian 
Higgs strings always reconnect, even at ultrarelativistic 
speeds (P = 1); unlike for other types of defects [231 135] . 
and against naive expectations, the only way in which 
strings can pass through each other appears to be by an 
even number of reconnections. 
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